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Abstract. We study the dynamics of quantum correlations of two coupled double 
quantum dots containing two excess electrons. The dissipation is included through 
the contact with an oscillator bath. We solve the Redfield master equation in 
order to determine the dynamics of the quantum discord and the entanglement of 
formation. Based on our results, we find that the quantum discord is more resistant 
to dissipation than the entanglement of formation for such a system. We observe that 
this characteristic is related to whether the oscillator bath is common to both qubits 
or not and to the form of the interaction Hamiltonian. Moreover, our results show that 
the quantum discord might be finite even for higher temperatures in the asymptotic 
limit. 
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1. Introduction 

Among the "unusual" manifestations observed in the quantum world, entanglement 
is undoubtedly one of the most interesting ones. Entanglement corresponds to global 
states of two or more quantum systems that cannot be separated into direct product 
states of individual subsystems. This characteristic yields correlations between quantum 
systems that cannot be found in any classical system pp. Moreover, entanglement 
is a very important ingredient for quantum computer's architecture [2] and quantum 
communication [1]. It is well-known that some tasks executed by a quantum computer 
can be performed exponentially faster than by the existing computers. These potential 
applications of a quantum computer have motivated a large number of experimentalists 
and different approaches have been used to build entangled states in laboratories 
[3]. Despite the successful experimental achievements, there are many difficulties to 
overcome before a functional quantum computer becomes a reality. One of the most 
trivial difficulties is the fact that the usual quantum systems candidates for qubits (see 
below) are intrinsically open to their environments and therefore suffer their effects 
by losing coherence. Thus, pure quantum states become mixed states under the 
environment's influence. Fortunately, it has been demonstrated that computers based 
on mixed states can also be used to solve certain tasks more efficiently than classical 
computers jl], although they are less powerful than the computation using pure states. 
The reason for such a performance is attributed to correlations not presented in classical 
systems. These quantum correlations (QCs) can be identified through a quantity 
called quantum discord [5]. Recently, the interest on this subject has received great 
attention due to the possibility of achieving quantum computation without entangled 
states [6]. As a consequence, many different aspects of quantum discord have been 
discussed, e.g., robustness to quantum sudden death [7J[8], relation to the speed-up in 
deterministic quantum computation with one qubit [9], the interplay between quantum 
phase transitions and quantum discord [10] . and the condition to obtain completely 
positive maps [TTj . 

There are many different systems proposed as candidates for qubits. Among them, 
double quantum dots (DQDs) (see Fig. [I]) are very interesting due to the easy integration 
with existing electronics and the advantage of scalability (T2j. Moreover, arrays of 
quantum dots can be employed to perform logical operations, which can be used to 
implement an universal quantum computation based on quantum [13] and classical 
effects [H]. In this work, we investigate the quantum correlation dynamics of two DQDs 
including the basic elements to simulate a realistic situation, e.g., finite temperatures, 
interaction between qubits, and detuning. We also compare the effects of dissipation on 
the behavior of the quantum discord (QD) with that of the entanglement of formation 
(EoF) [15] . Our results show the longer duration of QD against the EoF. To explain 
such a result, we analyze how the environment acts on the system by comparing the 
case where the bath is coupled to both DQDs, which is the more realistic case, to the 
case where each DQD is coupled to its own bath. Based on this comparison, we verify 
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Figure 1. (Color online) Schematic representation of two DQDs with one electron 
localized at the left (right) side of the top (bottom) DQD. See Ref. [T!5] for further 
details. 



that the common character of the bath might be very important to prolong the QD. 
In fact, this character of the environment brings different properties to the dissipative 
dynamics. Contrary to the case of independent environments, by which the coherences 
are completely lost, we can find regions of the state space where the coherences are 
preserved [16]. We observe that this robustness can be transferred to the quantum 
discord. 

The present paper is organized as follows. In Sect. II we describe the basic concepts 
of quantum discord. In Sect. Ill we present the Hamiltonian of the system and the master 
equation which accounts for the interaction between the system and the environment. 
In Sect. IV, we discuss the obtained results for the quantum correlations of two coupled 
DQDs. Finally, we summarize our work in Sect. V. 

2. Quantum Discord 

Quantum discord has been proposed as a quantity that captures the quantum 
correlations between two subsystems [5]. To evaluate such QCs, one needs to subtract 
the classical correlations from the total correlation [5j [17]. The total correlation of a 
bipartite quantum systems, A and B, is calculated by the quantum mutual information 



where S(A\B) = S(AB) - S(B) and S(X) = -TrpHogX} is the von Neumann 
entropy of the density matrix X. Here, we adopt the following definitions: AB = pab, 
A = pa = Tr B {pAB}, and B = p B = Tta{pab}, where pab is the bipartite density 
matrix. 

To calculate the classical correlations we observe that the projective measurements 
on a subsystem remove all nonclassical correlations between the parts, i.e., after a 



I(A : B) = S(A) — S(A\B), 



(1) 
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measurement on a particular subsystem, all QCs are destroyed. Thus, we define a 

quantity that evaluates the mutual information after a measurement on one of the 
subsystems 

l(A:B) {uf} = S(A)-S(A\{Uf}), (2) 

where S^AKII^}) is the system conditional entropy after the measurement and {nf } 
defines a complete set of one-dimensional projectors, and the different outcomes of 
this measurement are accounted for by j. To quantify the quantum correlations, 
since X(A : S)r n Bi depends on the projector basis {n^}, we take the maximum of 
T(A : B)mB\ taking into account all possible projectors. Thus, we can define the 
following quantity, 

J(A : B) = S(A) - win {uf} S(A\{Uf }), (3) 

that gives a measure of the total classical correlations between two subsystems p2] . 
Therefore, the QD can be written as: 

5(A: B)=1(A: B)-J(A: B). (4) 

It is straightforward to observe that if all information can be obtained locally by B, this 
subsystem has QD equals to zero. This implies that a measurement on B does not alter 
the state of A. Otherwise, if just part of the information can be obtained locally, B is 
quantum correlated to A. The QD thus gives the total information that is not locally 
accessible and has been accepted as a measure of the quantum correlation. Besides, the 
QD is possibly finite even for separable quantum states and can be utilized as a new 
resource for quantum computation [U [6] . Also, when only pure states are considered, 
QD and the entanglement of formation [15] are indistinguishable. 

3. Theoretical Model 

3.1. Hamiltonian of the System 

The system considered in this work is composed of two double quantum dots (see Fig. 1), 
where each DQD has an excess electron localized in either the left \L) or the right dot 
\R). This system can be modelled by the following pseudospin Hamiltonian 

H s = ^{a^+a^)+Ja^®af\ (5) 

where the first term describes the tunnelling coupling energy A between the qubits, 
and the last term takes into account the effects of the Coulomb interaction between 
the two electrons, which favors anti-parallel configurations \L, R) and \R,L) over the 
parallel ones \L,L) and \R,R). We adopt the convention \L) = | l) and \R) = 1 1)- For 
simplicity, we set the energy offset of the single-qubit states equals to zero. 
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3.2. Master Equation 

In order to analyze the two DQDs dynamics in an open quantum system, we suppose 
that both qubits are coupled to a bath of harmonic oscillators (phonons). The total 
Hamiltonian that computes the environment perturbation is given by 

H = H s +(a^+a^)C + H B , (6) 

where C = B + with B = J^kd^k, H B = J2k UJ k a l : a ky u k is the frequency of the 
k th normal mode of the bath , and h — 1. Here a& is the usual annihilation operator of 
this mode and gk is its complex coupling constant with the DQDs. It is important to 
note that the interaction Hamiltonian acts in a correlated way with both qubits, thereby 
representing a common bath. To numerically determine the evolution of the two DQDs 
reduced density matrix, we employ the Redfield master equation [20] 
dp At) 
dt 

where Hi(t) is the interaction Hamiltonian in the interaction picture, 

Hj(t) = UlWl(t) { (a« + C} U B (t)U s (t), (8) 

where U B (t) = exp (—iH B t) and Us(t) = exp (—iH$t). Explicitly, the unitary evolution 
Us(t) is given by 

U s (t) = - [cos (Qt) + cos (Jt)} I 

+ llcosm- cos (Jt)}a^aW 



- j dt'Trg {[Hj(t), [Hj(t'), PBPi(t)]}} , (7) 



sin ( Jt) — ^- sin (Qt) 
sin ( Jt) + — sin (Qt) 



y v 



0) 



where Q = a/J 2 + 4A 2 . Here, we suppose that the oscillator bath density matrix p B is 
initially decoupled from the system, 

p B = ^ex P (-(3H B ), (10) 

where Z is the partition function Z = Tr B [exp(— (3H B )], (3 = l/k B T, k B is 
Boltzmann constant, and T is the absolute temperature of the environment. Defining 
U B (t)CU B (t) = C(t) and A(t) = Ul(t) (ajp + a? ] ^j U s (t), we can write the interaction 
Hamiltonian in the interaction picture as follows: 

H I (t) = A(t)C(t). (11) 
Thus, substituting Eq. (ITUl) and Eq. (TTTj) in the master equation, Eq. (JT]), we obtain 

^= J*dt"D(t,t)[A(t), Pl (t)A(t')] 

+ f dt'V*(t,t')[A(t') Pl (t),A(t)}, (12) 
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where V{t, t') = T l (t- t') + T 2 (t - t') with 

T 1 (t-t') = Tr B {B(t)p B B\t')} 

= ^2 \9k\ 2 n k exp[-iu k (t - t')], (13) 
k 

and 

T 2 (t -t') = Tr B {B\t)~p B B(t')} 

= \9k\ 2 (n k + 1) exp[iu k (t - t% (14) 

k 

with B(t) = U B (t)BU B (t), p B = U B (t)p B U B (t), and n k is the average occupation 
number of mode k: 

n k = r. (15) 

exp(/3u k ) - 1 

Defining, as usual, the spectral function 

J(u) = ^2\g k \ 2 5(u - u k ), (16) 

we can replace the summations above by integrals 



CO 



^i(^) = / du J(u)n(u) exp(— iut), (17) 

and 

T 2 {t)= I dujJ(u))exp(icut)[n(u)) + l]. (U 



If we assume an ohmic spectral density to the reservoir, J(oo) = rju exp(—u/uj c ), where 
uj c is a cutoff frequency and rj is the damping constant, we can explicitly evaluate the 
integrals above which yield: 

P(M')- ^ 



[1 + iu c {t - t')] 2 

+ pRe {¥ lst \l + l/(p Ue ) - i(t - t')/(3)}, (19) 
where \l/( lst )(-) is the first polygamma function. 

4. Quantum Correlation Dynamics of DQDs 

To calculate the DQDs quantum correlation dynamics we numerically solve Eq. (j4]) using 
the density matrix dynamics computed by Eq. (J7J) as input. In our simulations we fixed 
A = nh/2r (A w 10/ieV) and the cutoff frequency u c = 200/r j2TJ, where r = 10" 10 s. 
These parameters correspond to typical experimental values in GaAs/AlGaAs lateral 
DQDs [19]. Our expressions are given as a function of the weight of the Coulomb 
interaction between the two electrons, J, and the damping constant i]. First, we analyze 
the unitary dynamics in Fig. [2] supposing that the two DQDs are initially uncorrelated, 
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Figure 2. (Color online) Quantum correlations (QD and EoF) for the unitary 
evolution ruled by H$. We consider different values of J/ A and we fix the coupling 
constant rj = 0. In (a), (b), and (c) we have J = A/8, J — A/4, and J — A/2 
respectively. This is what we define as weak coupling regime. For the strong coupling 
regime we have (d), (e), and (f) with J = 4A, J = 8A, and J = 16A respectively. 



Pi(fy = I Ti)(T4 I; an d assuming that the coupling with the environment is zero, i.e., 
7] = 0. Since we suppose an initial pure state, the EoF and the QD dynamics coincide. 
In Fig. [2] we observe that similar dynamics can be reached for different Coulomb 
coupling values J. Thus, one can find a entanglement dynamics in the weak coupling 
regime (J << A) that corresponds to a similar dynamics in the strong-coupling regime 
(J >> A), and vice versa. To determine these equivalent entanglement dynamics, we 
examine the period that the initial state | ti) reaches a maximum entangled state t max 
for J << A and J » A. For this purpose, we estimate these periods analyzing the 
coherences of Us(t)p(0)Ul(t). In the weak coupling regime J << A we obtain 

tmax « -J- (| + ni^j , (20) 

where n £ N and J< denotes that J « A. For the strong coupling regime J » A 
(</>), the period that the initial state reaches a maximum entangled state depends on 
A and is given by 

tmax - -^[^+n-) . (21) 
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Figure 3. (Color online) Two-qubit quantum correlations (QD and EoF) under a 
dissipative dynamics for J = 4A and T = 0.1K in (a), T = 0.5K in (b), and T = 2K 
in (c). The blue (dotted) is the QD dynamics and the red (solid) line is the EoF 
dynamics. Here we use 77 = 1/600. 



Based on Eqs. ( 1201) and ( l2TT) we find that the time to achieve a maximum entangled 
state is approximately the same in both regimes when J > J < = 2A 2 . This behavior can 
be noted by comparing Fig. [2]-a (J < = A/8) with Fig. [2]-f (J> = 16A) and Fig. [2]-b 
( J< = A/4) with Fig. He (J> = 8A). 

In Fig. [31 assuming that p/(0) = 1 1^) (t4- 1 j we present the numerical results for the 
QD and the EoF dynamics within the strong coupling regime (J = 4A), for different 
temperatures (T = 0.1K, T = 0.5K, and T = 2K) and for finite damping constant 
(77 = 1/600). The EoF for two qubits can be expressed in terms of the concurrence [22] . 
EoF = -/(C)log 2 (C) - (1 - /(C))log 2 (l - f(C)), where f{C) = (1 + s/T^)/2. 
On the other hand, the concurrence C is defined as the maximum between zero and 
Ai — A 2 — A 3 — A4, where Aj, for i = 1,2,3,4, is the square root of the eigenvalues 
of the matrix pa y ® a y p*a y ® a y , with Ai being the largest one among them, and p* 
is the complex conjugate of p. In Fig. [3] one notices that the QCs (QD and EoF) 
oscillate rapidly, which resembles the behavior observed when there is no coupling to 
the environment (see Fig. [2] (d)). However, when the coupling to the bath is included, 
the QCs decay as a function of time due to the loss of correlations between the qubits 
and the environment. Moreover, we verify that QD is more robust than EoF for higher 
temperatures and for higher values of the coupling to the environment rj (see Fig. HJ), 
and remarkably reaches a constant value in the asymptotic limit, when t — > 00. In Fig. [5] 
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Figure 4. (Color online) Two-qubit quantum correlations (QD and EoF) under a 
dissipative dynamics for J = 4A and T = 0.1K in (a), T = 0.5K in (b), and T = 2K 
in (c). The blue (dotted) is the QD dynamics and the red (solid) line is the EoF 
dynamics. Here we use 77 = 1/200. 

and Fig. [6l the EoF and QD dynamics are analyzed for different temperatures and bath 
couplings r) within the weak regime J = A/4. The QCs for 77 = 1/600 presented in 
Fig. [5] have a oscillation period that corresponds to the unitary evolution (Fig. 2(b)). 
However, when the coupling to the environment is more relevant rj = 1/200, we observe 
a short oscillation period (see Fig. [6]). Again the QCs decay as a function of time due to 
the interaction with the environment in both Figs. [5] and [6] By comparing the results of 
the weak coupling regime to those in the strong coupling regime, we notice that the QD 
is also more robust than EoF for higher temperatures and the QD tends to a finite value 
in this case. This behavior is related to the convergence of the system when t — > 00 to 
quantum subspaces defined by the eigenvectors of + ai 2 \ i.e., the system operator 
that is coupled to the environment [23]. In such a case, the resultant system density 
matrix can be written as 

p(t -+ 00) « PilttXttl +P 2 |*±)(*±| + Ps\U)(U\, (22) 

where |\l/±) = a\ t-i-)+/?| IT); an d Pi, P2, and are the weights of each eigenvector in 
the asymptotic limit. Because cxi 1 ^ + ai 2 ^ is degenerated, the density matrix coherence, 
whose modulus is given by p2\ct(3\, is sustained. This fact induces a finite quantum 
correlation that is preserved by the two DQDs. On the other hand, the equilibrium 
density matrix describes a disentangled state since p2\ct/3\ < ^P\Pz [21]. To verify that 
this is the cause of the robustness of the two DQDs quantum discord, we compare two 
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Figure 5. (Color online) Two-qubit quantum correlations (QD and EoF) under a 
dissipative dynamics for J = A/4 and T = 0.1 K in (a), T = 0.5 K in (b), and T = 2 
K in (c). The blue (dotted) is the QD dynamics and the red (solid) line is the EoF 
dynamics. Here we use rj = 1/600. 
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Figure 6. (Color online) Two-qubit quantum correlations (QD and EoF) under a 
dissipative dynamics for J = A/4 and T = 0.1 K in (a), T = 0.5 K in (b), and T = 2 
K in (c). The blue (dotted) is the QD dynamics and the red (solid) line is the EoF 
dynamics. Here we use rj = 1/200 




Figure 7. (Color online) EoF (a) and QD (b) as a function of time, considering J = 
for independent (dotted curve) and common (solid curve) baths with T — 0.1 K and 
r/ = 1/600. The insets show the results in logarithmic scale. 



distinct situations. In the first situation, each DQD is coupled to its own heat bath, 
which are completely uncorrelated. The second situation has a common environment, 
i.e., the same bath is coupled to both DQDs and the interaction Hamiltonian is the one 
already considered + C. 

For independent environments, the interaction Hamiltonian in Eq. ([6]) is slightly 
different and is given by <Jz£^ + a z C^- Therefore, the master equation's 
correlation functions, Eq. ([13]) and Eq. (fl4l) . are written as T\{t — t') = 
^B{Et k/=1 B^(t)psB^\t')} andT 2 (t-t') = Tr^^i B^\t)p B B^(t')}, 
respectively. In this case, because the qubit's environments are totally uncorrelated, 
we have Ti(t — t') = T 2 (t — t') = if k ^ k' . We also suppose a simple case where the 
Coulomb coupling is null J = 0. Initially, the two DQDs are considered in a maximal 
quantum correlated state given by pi(0) = |$)($|, where |<3>) = (HJ,) + \l^))/\/2. 

By examining Fig. [TJ we observe that for independent environments the EoF and 
the QD exponentially tends to zero. On the other hand, for common environments, the 
QD has an exponential decay that tends to a finite value due to the specific form of the 
coupling between the system and the environment. The fact that the environment acts 
correlated with both DQDs makes a more robust QD. To investigate this robustness, in 
Fig. [8] we compare the EoF and QD dynamics considering only the common bath case 
and observe that QD is more resistant to the environment's disturbance. Moreover, the 
inset of Fig. |8]shows only the exponential decay of EoF and the QD, i.e., the EoF curve 
and the QD curve minus the long time saturation value. In this case, we notice that the 
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Figure 8. (Color online) Logarithmic scale of EoF (solid curve) and QD (dotted 
curve) as a function of time, considering J = and a common bath with T = 0.1K 
and T) = 1/600. In the inset we plot the EoF curve and the QD curve minus the long 
time saturation value to extract only its exponential behavior. 



exponential behavior of QD decays faster than EoF. 

Based on our results, we conclude that the quantum correlation dynamics of 
the two DQDs can be summarized as follows: the Coulomb interaction introduces 
the correlations between both DQDs and the environment acts in a way to preserve 
these correlations for a longer time. Despite the fact that the entanglement decays 
exponentially to zero, the QD is sustained depending on the way the DQDs are coupled 
to the environment. 

5. SUMMARY 

We have numerically solved the Redfield master equation to study the quantum 
correlation dynamics of two DQDs at finite temperatures including the interaction 
between the qubits and detuning. We have verified that dissipation in this system 
has an impact on the EoF which is stronger than on the QD, thereby causing a shorter 
duration of the EoF. We have shown that the QD is not completely destroyed through the 
interaction with the environment and that it remains stable even at finite temperatures. 
To explain these facts, we have explored the characteristics of the bath, i.e., whether 
both qubits are coupled simultaneously to a single bath or each qubit is coupled to its 
own bath. The results have shown that even though the coupling between the two-qubits 
is zero, the quantum correlations survive for longer periods in the case of a common 
bath. This behavior results from the particular form of the interaction Hamiltonian and 
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there are systems for which the robustness of QD may not occur and each different case 
must be analyzed carefully. 

As a final comment, it should be desirable to study the same effects in the context 
of the dynamics of a couple of Brownian particles in a common bath which has recently 
been analyzed in (25] and where it has been shown that the reservoir mediates an 
effective coupling between the particles. The alternative (and more general) coupling to 
the environment used in this problem could also be applied to the present model aiming 
at a possible maintenance of the quantum correlations between the qubits for a much 
longer time interval. 
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